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ABSTRACT 

We use a quadratic estimator with KL-compression to calculate the angular power 
spectrum of a volume-limited Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) 
galaxy sample out to I = 200. We also determine the angular power spectrum of 
selected subsamples with photometric redshifts z < 0.3 and 0.3 < z < 0.4 to examine 
the possible evolution of the angular power spectrum, as well as early-type and late- 
type galaxy subsamples to examine the relative linear bias. In addition, we calculate 
the angular power spectrum of the SDSS DR7 main galaxy sample in a ~ 53.7 square 
degree area out to I = 1600 to determine the SDSS DR7 angular power spectrum to 
high multipoles. We perform a x 2 fit to compare the resulting angular power spectra 
to theoretical nonlinear angular power spectra to extract cosmological parameters and 
the linear bias. We find the best-fit cosmological parameters of O m = 0.267 ± 0.038 
and n b = 0.045 ± 0.012. We find an overall linear bias of b = 1.075 ± 0.056, an early- 
type bias of b e = 1.727 ± 0.065, and a late-type bias of 6/ = 1.256 ± 0.051. Finally, we 
present evidence of a selective misclassification of late-type galaxies as stars by the 
SDSS photometric data reduction pipeline in areas of high stellar density (e.g., at low 
Galactic latitudes). 
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1 INTRODUCTION 

The amount of normal, baryonic matter, dark matter, and 
dark energy in the Universe controls the formation and evo- 
lution of large scale structure. Hence, by observing the dis- 
tribution of the large scale structure, as traced by the galaxy 
density, we can determine the densities of these components 
of the Universe. The distribution of galaxies can be statis- 
tically qu antified by the three dimensional power spectrum 
(see, e.o.. |Percival et aZll2007l l. Spectroscopic surveys, how- 
ever are generally limited by either sample size, cosmic vol- 
ume, or both, as obtaining spectroscopic redshifts for a large 
sample of faint galaxies can be expensive and difficult. In- 
stead, we can utilize the large photometric surveys of galax- 
ies to more precisely measure the angular power spectrum 
(|Tegmark et qZj|2002h . and compare these measurements to 
a theoretical 3D power spectrum by using the photometric 
redshift distribution of the galaxy sample. 

From these comparisons, we can extract constraints 
on any p arameters that effe ct the power spectrum, includ- 
ing bias (iBlake et a~L 20041 '). cosmological mat ter density, 
flm dHuterer et al\ 1200 li b barvon density, flh dFrith et al\ 
120051 ). or even neutrino masses (|de Putter et al. 20121 ). 
In addition to the parameter space explored by galaxy 



angular power spectra, the natural relationship between 
the angular power spectrum and the 3D power spec- 
trum can be used to recon struct the 3D power spectrum 
(e.g., iDodelson etd\ |2002h . Due to these useful proper- 
ties, galaxy angular power spectra ha ve been estimated for 
many surveys, incl uding the SPS S dTegmark et all 12002;; 
[ Thomas et fflZl20ld:lHo" et aZj|2012h . EDS GC dHuterer et all 
I2001D , NVSS lBlake et aZj2004h 2MASS dFrith et al.ll2005t ). 
and others (e.g., Baugh fc Efstathioulll994 ). 

In a previous paper dHaves et al\ 120121 : hereafter 
HBR1 2), we detailed our implementation of the lBond et all 
( 1998) quadratic estimator to calculate the angular power 
spectrum of the SDSS DR7 using HEAhPix. With this tech- 
nique, we used our calculated angular power spectra to- 
gether with a projected linear 3D power spectrum to re- 
trieve constraints on cosmological matter density and bias. 
The cosmological parameters of greatest interest, however, 
impact the angular power spectrum at multipoles that cor- 
respond to small scales, that is in the nonlinear regime. To 
more tightly constrain cosmological parameters, it is neces- 
sary to fit estimated angular power spectra to theoretical, 
nonlinear 3D matter power spectra, for example, by using 
the Code for Anisotro pies in the Micro wave Background, 
also known as CAMB (|Lewis et aZ.ll2000l V 
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In this paper, therefore, we build on this previous work 
to present new angular power spectrum measurements. How- 
ever, we have chosen first to calculate angular power spectra 
of galaxies from the SPSS DR7 using a volume-limited sam- 
ple created bv lRoss et oil (|2010h . This allows us to make sub- 
samples that can be directly compared to each other, such as 
high and low redshift samples to investigate the evolution of 
the angular power spectrum, as well as galaxy-type defined 
samples that allow for a measurement of the relative lin- 
ear bias between early- and late-type galaxies. This current 
paper first begins with an overview of our data, quadratic 
estimator, and x 2 fitting process in Section [5] Next, we 
use CAMB to determine theoretical angular power spectra 
to compare with the estimated SDSS DR7 angular power 
spectra and detail the effects this has on the \ 2 fits to ex- 
tract cosmological parameters. In Section [4] we determine 
the galaxy angular power spectrum using a volume-limited 
sample of SDSS DR7 galaxies that is further divided into 
two redshift slices to quantify the redshift evolution of our 
quantified cosmological parameters and the linear bias. We 
split this volume-limited sample into early- and late-types in 
Section[5] and perform the same analysis to examine the rel- 
ative bias between galaxy types, as well as combine the type 
cuts and redshift slices to look at the evolution of the rela- 
tive bias. Finally, we examine the angular power spectrum 
of SDSS DR7 galaxies at high multipoles, by constraining 
our analysis to a subset of the full DR7 and applying our 
quadratic estimator to a higher resolution map in Section [6] 
We discuss the results in Section [7] and conclude in Section 

m 



2 OVERVIEW 

The basic data, systematic cuts, and angular power spec- 
trum estimation technique are fully detailed in HBR12. Here 
we provide a concise overview to explain the basics of our 
analysis. 



2.1 SDSS DR7 Data 

We have chosen a contiguous area of the SDSS DR7 North 
Galactic Cap ellipsoid to examine, using SDSS stripes 9-37. 
Our SDSS DR7 data consists of all galaxies in this area with 
apparent r-band magnitudes from 18-21 with photometric 
redshifts and photometric redshift errors. We have created 
subsamples of this data for our analysis, a volume-limited 
sample and associated subsamples described and used in 
Sections [4] and [5] and a selected area of the main galaxy 
sample around the North Galactic Pole described in Section 

El 

These samples are pixelated with HEALPix and we 
mask pixels that fully or partially lie outside of the SDSS 
DR7 footprint, have average seeing greater than 1.5 arcsec- 
onds, have average reddening greater than 0.2 magnitudes, 
and have areas of poor observing quality. Furthermore, any 
pixel that has less than 75% usable area after masking is 
also rejected. After pixelating with HEALPix and masking, 
we create a vector of pixelized overdensities: 



(1) 



where Gi is the number of galaxies in pixel i, G is the survey 
average galaxy density, and fii is the pixel area. 

As the SDSS is not a full sky survey, we are prevented 
from determining each multipole moment individually. Mul- 
tipole resolution is determined by A£ ~ 180° /<j>, where <j> is 
the smallest angular dimension of the survey (|Peebleslll980l l. 
By using the full area of the SDSS, we can achieve a mul- 
tipole resolution of A£ « 5, and thus we group multipole 
moments together into bands. Since we cannot determine 
the power of individual multipole moments, we have as- 
sumed that multipole moments in a band have the same 
value, which also simplifies the computational complexity of 
this quadratic estimation method. Since we do not expect 
large variations between adjacent multipoles in the galaxy 
angular power spectrum, as evidenced by the smooth the- 
oretical predictions, this assumption is both justified and 
convenient. 



2.2 Quadratic Estimation 

For our implementation of the quadratic estimation tech- 
nique, we first construct a covariance matrix from our pix- 
elixed overdensities, and we want to compare this to a model 
covariance matrix composed of a signal matrix plus a noise 
matrix: 



(xixj) = Sij + N K 



(2) 



where the signal matrix describes the covariance matrix we 
expect based on an initial angular power spectrum and the 
geometry of our survey area, and the noise matrix is modeled 
as a Gaussian random process: 
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where the Pe are the Legendre polynomials, fty is the angle 
between pixels i and j, t is the beam width, and Oi is the 
rms noise in pixel i. Using the assumption that the Ci in a 
band are equal, we define the Pi, matrices for each band: 
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(5) 



Our goal is to determine the Ct that produce a model 
covariance matrix that matches the covariance matrix con- 
structed from our pixelized data. We h ave im plemented 
the quadrati c estimator of iBond et al\ (|l998l) with KL - 
compression ijVoeelev fc Szal av 1996; Tcg mark et al.lll997f ): 



SC b 



:( F ' 



Tr[(x , x' T -N')(C'- 1 P;,C'- 1 )] 



where, the Fisher information matrix F is defined as: 



F W = ~ Tr ( C'-MC'^PI/ 



(6) 



(7) 



We u se F 1//2 for this estimator as suggested by 
iTegmarkl (1 19981 ) to provide uncorrelated bandpower error 
bars, given by o\> = \J (F~ 1 )bb, and well behaved win- 
dow functions. This computationally demanding process can 
be iterated and converges quickly given a reasonable input 
power spectrum. 
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After determining the angular power spectrum, we can 
recombine adjacent bandpowers together to produce higher 
signal-to-noise measurements at the cost of reducing multi- 
pole resolution. We determine the midpoint of each band as 
the point where half the power in the band is below the mid- 
point and half above, and the endpoints of each band where 
the power falls to e _1//2 of the peak power. This step is useful 
for plotting the results, though we use the non-recombined 
angular power spectra in our \ 2 fits. 

2.3 x 2 Fitting 

After determining the angular power spectra of our data, 
we want to compare our observed results to the results 
predicted by a given theory, or alternatively to determine 
the theoretical parameters that best fit our data. We do 
this by constructing a theoretical angular power spectrum 
from a prior 3D power spectrum that is dependent on cos- 
mology. We then project this 3D power spect rum down to 
two dimensions usi ng Limber's approximation IjLimber 1953; 
ICrocce et fflHl2010t ) and the redshift distribution of the sam- 
ple under consideration: 

Ce " wtt) J * {z)D {z)p( ^r ) ^) b dz (8) 



where D(z) is the growth function (|Carroll et al\\W9$ ) . b is 
the bias, and r and ~g are the comoving distance and number 
density respectively. 

With the theoretical angular power spectrum, we then 
can compare with the observed angular power spectrum 
through x 2 fitting: 

X 2 K) = J^OnCi - lnC 6 T ) c »FwC b , (lnC v - InCj) (10) 

bb' 

By minimizing this \ 2 , we can determine the 3D power spec- 
trum that best fits the data and thus can constrain the pa- 
rameters that produce that power spectrum. 



3 FITTING TO THEORY 
3.1 Nonlinear Power Spectrum 

To improve the precision of our measured constraints on 
cosmological parameters over our previous fits, we imple- 
mented a theoretical angular power spectrum calculation us- 
ing nonlinear 3D power spectra obtained by using the Code 
for Anisotropies in the Microwave Background (CAMB: 
iLewis et al\\200dl ). Among its many functions, CAMB is 
capable of produ cing nonlinear 3D matter power spectra 
using HALOFIT jSmith et al.||2003h . HALOFIT generates 
nonlinear matter power spe ctra using the halo model to 
descr ibe galaxy correlations dSeliak 2000; P eacock fc Smith! 
2000). The halo model suggests that the collisionless dark 
matter forms dark matter haloes within which baryons col- 
lapse to form some number of galaxies. Thus, the matter 
power spectrum is separated into two regimes: the corre- 
lation of galaxies is determined by the correlation between 
different dark matter haloes on large scales, and the cor- 
relation between galaxies within the same haloes on small 



scales. HALOFIT has determined empirical fitting functions 
for the power spectra of these two regimes by matching N- 
body simulations over a range of cosmological parameters. 
The sum of the quasi-linear large scale term and the non- 
linear small-scale term produces an overall nonlinear matter 
power spectrum. 

By running CAMB and HALOFIT with high accuracy 
options enabled, we produce no nlinear matter powe r spec- 
tra with an accuracy of ~ 0.2% (|Howlett et qZj|2012h . How- 
ever, this accuracy does assume a number of factors, includ- 
ing: the ionization history of the Universe, the Hubble pa- 
rameter, the dark energy equation of state parameter, and 
the initial power spectrum. Smaller issues in these calcu- 
lations include a number of other parameters such as the 
primordial Helium fraction and effective number of neutrino 
species. Where appropriate, we assume a fiat cosmology and 
WMAP 7-year best fit results and leave other parameters 
at the CAMB default settings, which correspond to the cur- 
rent best measurements of these parameters. We generate 50 
matter power spectra from z — to z = 1 with Az = 0.02 
for each set of Q m , fit, and bias b that are fit, and proceed 
in the calculation from Equation [8] as before. 

3.2 Nonlinear Fits 

Now that this calculation can be extended to nonlinear 
scales, we can determine cosmological parameters more pre- 
cisely than in HBR12. In addition to the bias and Q m , we 
also chose to fit the baryon density, Qf We also investigated 
fitting the spectral index, n s , but found that generally the 
spectral index altered the angular power spectra similarly 
to fi m , thus causing a degeneracy in the parameter fits. As 
a result, we chose not to fit this parameter. 

Before fitting, we must determine the maximum multi- 
pole to which we can fit the estimated and theoretical an- 
gular power spectra. In some cases, we are limited by the 
signal-to-noise of the data, as we can't fit into the noise 
dominated region of our observed angular power spectra. 
But in most samples, the signal-to-noise is sufficient for the 
full range of £ in our estimation and we are instead limited 
by the pixelization scale. The pixelization process causes a 
loss of information on scales near the pixel size, but this 
is not a sharply defined boundary. Instead, the pixeliza- 
tion discards steadily more information as you approach the 
pixel size. The power lost due to pixelization is given by 
the pixel window functions, we, calculated and supplied by 
HEALPbJE which is demonstrated in Figure [T] for resolution 
64 (i.e. HEALPix NSIDE parameter 64), the resolution of 
our volume-limited sample and subsample results. We see 
that by £ ~ 175, the pixelization supresses 50% of the power 
in the angular power spectrum, even though the equivalent 
linear scale of a pixel is £ ~ 200. The pixel window func- 
tions relate the observed pixelated angular power spectrum 
C^ tx and the underlying unpixelized angular power spectrum 
C e " by: 

C% =w e C e y (ll) 

This has the effect of strongly supressing the theoretical an- 
gular power spectrum at high £, but also reducing power at 

1 http:/ /healpix. jpl.nasa.gov/html/intronodel4. htm 
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Multipole I 



Figure 1. The pixel window function at HEALPix resolution 
64 in black, demonstrating the percentage of power lost at each 
multipole for any angular power spectrum calculated at this pix- 
elization scale. In red, we integrate the pixel window function up 
to I to show the cumulative power lost to that multipole. We also 
show the maximum multipole we use in our \ 2 fits at I = 100, 
equivalent to twice the linear scale of the pixel. 



all scales. We have chosen to fit to a maximum multipole 
equivalent to twice the pixel scale for all nonlinear theory 
fits, where the bandpower value is reduced by roughly 20% 
by the pixel window function, but is still dominated by the 
unpixelized angular power spectrum. 

Finally, we must be aware of a known systematic 
|Tegmark et all < |2002h > in this quadratic estimation method 
that causes an excess of power at the small scale (i.e. high 
£) end of the angular power spectrum. Typically, this only 
affects the last few bandpowers in our estimation; and to 
correct for this effect, we use the quadratic estimator to cal- 
culate out to scales equivalent to ~ 80% of the pixel size, and 
discard the extra bandpowers. However, due to covariance 
between the bandpowers, some of this power is inevitably 
transferred onto small scales, complicating the exact calcu- 
lation. 



4 VOLUME-LIMITED SAMPLE 

For the entire SDSS DR7 galaxy sample, we are magnitude- 
limited, which biases the sample toward detecting intrin- 
sically brighter galaxies at higher redshifts. To improve 
on this, we construct a volume-limited sample out to red- 
shift z = 0.4 to examine a complete sample that is rela- 
tively free from M almquist bias. Th is volume-limited sam- 
ple was created bv lRoss et all (|2010t ). selecting over 3.2 mil- 
lion DR7 galaxies with de-reddened r-band apparent mag- 
nitudes m r < 21 to ensure galaxy completeness, r-band ab- 
solute magnitudes M r < —21.2, and photometric redshifts 



1.5 




0.1 0.2 0.3 0.4 0.5 



Multipole I 

Figure 2. The redshift distributions of the redshift samples with 
the z < 0.3 sample in black and 0.3 < z < 0.4 in red. These 
samples cover approximately equal cosmic volumes; but the high 
redshift sample is more peaked than the low redshift sample. 

z < 0.4. This absolute magnitude cut is more stringent 
than the SDSS detection limit to account for differences in 
k-corrections between early- and late-type galaxies. These 
data are also masked for seeing greater than 1.5 arcseconds, 
reddening greater than 0.2 magnitudes, and areas of poor 
image quality, following the method outlined in HBR12. 

4.1 Volume-Limited Data 

This volume-limited data set contains 3.2 million galaxies 
between redshift 0.0 < z < 0.4. We have additionally sep- 
arated this volume-limited sample into two mutually exclu- 
sive redshift slices of approximately equal cosmic volume 
from z < 0.3 and 0.3 < z < 0.4 to examine the possible 
evolution of the angular power spectrum with redshift and 
possible variation of the fit constrained cosmological param- 
eters. The redshift distributions of the two redshift slices is 
shown in Figure [2] 

4.2 Volume-Limited Angular Power Spectrum 

We have calculated the angular power spectra of our low 
redshift z < 0.3 and high redshift 0.3 < z < 0.4 sam- 
ples and show them in Figure [3] We can see that generally 
these two angular power spectra look roughly indistinguish- 
able, though perhaps there is some variation at multipoles 
100 < £ < 150. To examine this closer, we have taken the 
ratio of these angular power spectra, which is shown in black 
in Figure [4] We see that the ratio is fairly consistent with 
a value of one over the range of theoretical fitting, £ < 100, 
which implies that we cannot confidently detect any signifi- 
cant evolution at these scales. Beyond the range that we fit 
a theoretical angular power spectrum to, 100 < £ < 150, we 
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Angular Power Spectrum Ratio 
Ratio z<0.3 to 0.3<z<0.4 




100 
Multipole I 



Figure 3. The best fit theoretical angular power spectrum to 
the samples split up by redshift shell with the z < 0.3 sample 
in green and the 0.3 < z < 0.4 sample in red. Also apparent in 
the theoretical curves for the high redshift sample, plotted for 
£ < 100, are the wiggles from baryon acoustic oscillations. The 
BAOs are smoothed over in the low redshift sample due to the 
larger redshift range. 



do see a small deviation from unity, slightly outside the la 
error bars, which is not reflected in our fit values. However, 
beyond I = 150, we see that the values return to consistency 
with unity. 

Recall that although our samples are of galaxies with 
measured photometric redshifts z < 0.3 and 0.3 < z < 0.4 
and are therefore mutually exclusive, each of these galaxies 
has a photometric redshift error associated with them. We 
assume a Gaussian probability distribution function for each 
galaxy with mean equal to the measured photometric red- 
shift and standard deviation of the m easured photometric 
redshift error. As in Ro ss et ~al\ (|2010t ). we take 10 samples 
from the PDF of e ach galaxy, which are then weigh ted by lu- 
minosity function (|Montero-Dorta fc Pradal ( 2009)) and vol- 
ume constraints to account for the lower likelihood of both 
the galaxy being bright and residing in a smaller cosmic vol- 
ume (i.e. lower redshift). The galaxy redshift samples are 
normalized so that each galaxy contributes equally to the 
estimated "true" redshift distribution, which is shown in 
Figure [2] for the redshift subsamples of the volume-limited 
data set. Since these redshift distributions overlap, we ex- 
pect covariance between the measurements in each sample. 
Therefore, we expect some similarities in the features of the 
angular power spectra. 



4.3 Fitting Volume-Limited Results 

We have also produced fits out to a maximum of £ = 100 
(twice the pixel size) of the nonlinear theoretical angular 



Figure 4. The ratio of the angular power spectra of the redshift 
samples z < 0.3 and 0.3 < z < 0.4. There is a slight dip for 
100 < £ < 150 where the high redshift sample has more power, 
but generally the results are consistent with a ratio of one. 



power spectra to the volume-limited sample and subsamples. 
Due to the power associated with a high stellar density in 
the first bandpower of the late-type galaxy angular power 
spectra (see, e.g., Figure |S}, discussed in Section 17.11 we 
have not included the first bandpower in these fits. To be 
consistent, we have excluded the first bandpower from all fits 
derived from the volume-limited sample, though the effect 
on the non-late-type galaxy samples was small, changing 
the best-fit parameters by ~ 0.001. These fits are generally 
consistent with the results in Section 16.21 and the WMAP 
7- year results and show no strong evidence of evolution with 
redshift. These best fit theoretical spectra shown in Figure 
[3] Of int erest to note is that the baryon acoustic oscillations 
(BAOs) (|Eisenstein et ad 120051 ; ISeo et alihoi'j ) are clearly 
visible in the 0.3 < z < 0.4 sample due to its narrow redshift 
range, whereas the BAOs are smoothed over in the larger 
redshift ranges of the z < 0.3 and z < 0.4 samples. The 
best fit values of bias and cosmological parameters for these 
samples are given in Table [T] 



5 GALAXY MORPHOLOGY 

We also separate the volume-limited sample into early- and 
late-type galaxies, and, in addition to clearly showing the 
stronger clustering of early-type galaxies, we are able to ef- 
fectively determine the relative linear bias between these two 
morphological types in Section 15.11 These samples are cre- 
ated by separating early- and late-type galaxies based on the 
pztype parameter provided by the SDSS. The pztype pa- 
rameter is an estimate of the spectral type of the galaxy cal- 
cuated in the photometric redshift pipeline (|Abazaiian et all 
2009), and we classify galaxies with a pztype value of less 
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Redshift 


Type 


Bias 




Q b 


Full-z 
Full-z 
Full-z 


All 

Early 
Late 


1.545 ± 0.057 
1.727 ± 0.065 
1.256 ±0.051 


0.282 ±0.026 
0.275 ± 0.025 
0.262 ± 0.030 


0.041 ±0.020 
0.033 ± 0.019 
0.039 ±0.023 


Low-z 
Low-z 
Low-z 


All 

Early 
Late 


1.421 ± 0.054 
1.666 ± 0.064 
1.172 ± 0.050 


0.300 ± 0.036 
0.351 ±0.038 
0.279 ±0.041 


0.026 ±0.024 
0.034 ± 0.022 
0.032 ± 0.032 


High-z 
High-z 
High-z 


All 

Early 
Late 


1.515 ±0.057 
1.634 ±0.062 
1.363 ± 0.062 


0.234 ±0.019 
0.239 ±0.019 
0.244 ± 0.024 


0.036 ± 0.014 
0.033 ±0.014 
0.036 ± 0.017 



Table 1. The best fit biases and cosmological parameters for the 
volume-limited sample and subsamples. 



than 0.1 as earl y-type and those with pztype greater than 
0.1 as late-type |Ross et al\\20ldl ). 

In Section [521 we combine the redshift and galaxy type 
cuts together to produce four more subsamples that allow 
us to examine the possible redshift evolution of the angular 
power spectrum of early- and late-type galaxies. 




5.1 Galaxy Morphology Results 

As early-type galaxies are believed to preferentially form in 
high density environments, we wish to explore evolution in 
the angular power spectra of galaxies split by galaxy type. 
In our volume-limited sample, we can compare the angular 
power spectra of early-type galaxies, which tend to be larger 
and brighter, to late-type galaxies. This is important since 
by using a volume-limited sample we avoid Malmquist bias, 
which would result from selectively detecting only brighter 
galaxies in a magnitude-limited sample, which have more 
power at all scales. 

We also want to estimate the linear bias of early-type 
and late-type galaxies, and calculate the relative bias be- 
tween them. The results of our galaxy morphology angu- 
lar power spectra are given in Figure [5] Since the relative 
linear bias is roughly the ratio between these two angular 
power spectra, we can easily estimate the relative bias be- 
tween early- and late-type galaxies by examining the ratio 
in black in Figure [5] We can see that the relative bias is 
remarkably consistent across all but the very largest scales 
down to ~ 1 degree, and the bias of early-type galaxies is 
roughly 30 — 40% greater than the bias of late-type galaxies. 
The best fit values of bias and cosmological parameters for 
these samples are presented in Table [T] 



Figure 5. The angular power spectra of the early- and late-type 
galaxies in the volume-limited sample. 



Angular Power Spectrum Ratio 



Ratio Early to Late - z<0.4 
Ratio Early to Late - z<0.3 

Ratio Early to Late - 0.3<z<0.4 




50 100 
Multipole I 



5.2 Combining Redshift and Galaxy Morphology 
Cuts 

The results of our estimation of the angular power spectra 
of the different redshift slices for early- and late-type galax- 
ies are given in Figures [7] and [8] We see, as in Section I4.21 
no clear evidence of evolution between these two samples in 
different redshift ranges. Note that the high redshift late- 
type results are also signal-to-noise limited around t ~ 120, 
we therefore truncate that angular power spectrum. The ra- 
tio of angular power spectra of these samples are plotted 
along with the results for all types above in Figure [9l and 



Figure 6. The ratio of the angular power spectra of the early- 
to late-type galaxies in the volume-limited sample. 

are also generally consistent with a value of one suggesting 
no evidence of evolution. 

Finally, we want to look for the possible evolution 
of early- and late-type galaxies between our two redshift 
samples. By comparing the linear biases of the early- and 
late-type galaxies, we can also determine the relative lin- 
ear bias of these galaxy types. We find that the relative 
bias b e /bi — 1.375 ± 0.076 for our entire z < 0.4 volume- 
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Figure 7. The best fit theoretical angular power spectrum to 
the early-type samples split up by redshift shell with the z < 
0.3 sample in green and the 0.3 < z < 0.4 sample in red. Also 
apparent in the theoretical curves for the high redshift sample 
are the wiggles from baryon acoustic oscillations. The BAOs are 
smoothed over in the low redshift sample due to the larger redshift 
range. 



Figure 8. The best fit theoretical angular power spectrum to the 
late-type samples split up by redshift shell with the z < 0.3 sam- 
ple in green and the 0.3 < 2 < 0.4 sample in red. Also apparent in 
the theoretical curves for the high redshift sample are the wiggles 
from baryon acoustic oscillations. The BAOs are smoothed over 
in the low redshift sample due to the larger redshift range. 



limited sample, while b e /bi = 1.421 ± 0.083 for z < 0.3 and 
b e /bi = 1.200 ± 0.071 for 0.3 < z < 0.4. These are fairly 
consistent across all samples, and show no evidence of sig- 
nificant scale dependence on these large scales. 



6 LARGE MULTIPOLES 

We would like to determine the angular power spectrum 
to the highest multipoles (and thus smallest scales) allowed 
by the SDSS DR7 galaxy data, and to examine these small 
scales requires a high resolution pixelization. However, we 
are computationally limited by t he O(ni) sca ling depen- 
dence of the quadratic estimator l|Borrilll | l9 99) , where n p 
is the number of pixels in the map. By using large shared- 
memory supercomputers, this limits the analysis to maps 
with rip ~ 10 4 ; thus, in order to go to higher resolution 
and multipoles in the absence of better computational meth- 
ods, we are forced to constrain our analysis to smaller areas 
than the full SDSS DR7. Performing this calculation with 
a smaller area necessitates using larger bands as the band 
width is inversely proportional to the smallest linea r dimen- 
sion of the area under consideration l)PeeblesHl980r ). There- 
fore, obtaining the angular power spectrum at high multi- 
poles involves sacrificing band resolution and survey area, 
which means losing information about the angular power 
spectrum at all scales. However, this sacrifice is offset by be- 
ing able to constrain the angular power spectrum at smaller 
angular scales where cosmological parameters have a greater 



effect, ideally this will allow tighter constraints on those pa- 
rameters. 



6.1 Large Multipole Results 

To extend the angular power spectrum out to high multi- 
poles, we need high signal-to-noise and therefore we use the 
full magnitude- limited data of HBR12 instead of the volume- 
limited subset. This includes all galaxies with r-band mag- 
nitudes in the range 18-21 that have photometric redshifts 
and associated errors. 

When restricting our analysis to a smaller area, we are 
at least free to choose the particular area we examine. We 
select an area of low reddening near the North Galactic Pole, 
and limit our analysis to the ~ 53.7 square degree area corre- 
sponding to nested pixel number 162 at HEALPix resolution 
8. For the calculation, we use ~ 0.013 square degree pixels at 
HEALPix resolution 512, which allows us to compute bands 
of width 25^ out to £ ~ 1600, the scale roughly equivalent to 
the pixel size. This sample has a mean redshift of z = 0.251. 
We present the results from this measurement in Figure [101 



6.2 Fitting Large Multipole Sample 

Using the nonlinear power spectra produced by CAMB, we 
have fit the results of our large multipole sample for the 
parameters f2 m , fij,, and bias. We have fit out to a maximum 
of I — 800 or twice the pixel size to avoid signal lost by the 
pixelization. We find that the best fit Q m = 0.267 ± 0.038, 
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Angular Power Spectrum Ratio 

Ratio z<0.3 to 0.3<z<0.4 - Early 
Ratio z<0.3 to 0.3<z<0.4 - Late 





Multipole I 



Figure 9. The ratio of the angular power spectra of the redshift 
samples z < 0.3 and 0.3 < z < 0.4, separated by type. The 
late-type high redshift sample is signal-to-noise limited beyond 
£ = 120 and we truncate the angular power spectrum for that 
sample. We see a slightly more pronounced dip for 100 < £ < 150 
in the early-type galaxies compared to Figure|4] but we still can't 
strongly conclude that there is evidence of significant evolution 
in redshift. 



Q b = 0.045 ±0.012, and b - 
in Figure [TO] 



1.075 ± 0.056, and show this fit 



7 DISCUSSION 

7.1 Late-Type Large Scale Power 

In Figures [5] and [6l it is clear that the power in our first 
bandpower for late-type galaxies is unusually high, and this 
is consistent across all late-type samples. The smallest scale 
probed by this range of £ is over 30 degrees, a great deal 
larger than where we expect significant structure to exist. 
This suggests that despite our masking process correcting 
for seeing, reddening, and poor observing quality, there may 
be a systematic that preferentially affects late-type galax- 
ies on large scales. We have examined our late-type sample 
in more detail, and have found that the average density of 
late-type galaxies drops at either end (in lambda, the sur- 
vey longitude) of the SDSS observing footprint, nearer to 
the Galactic plane. While we already discussed in HBR12 
that the overall galaxy overdensity was consistent with zero 
at all Galactic latitudes, the late-type sample is more signif- 
icantly affected by this large scale systematic effect and this 
is apparent in the first bandpower of these angular power 
spectra. 

As this effect occurs closer to the Galactic plane, sys- 
tematics that could reduce late-type galaxy counts include 
stellar obscuration of background galaxies, variable star- 
galaxy separation efficiency, and an insufficiently strict red- 



Figure 10. In black, we show the angular power spectrum out 
to £ = 1600 of an ~ 53.7 square degree area of the SDSS DR7 
near the North Galactic Pole, at an eighth of the linear scale 
of the full area results. The best fit theoretical angular power 
spectrum to the large multipole, high resolution sample out to 
£ = 800, equivalent to twice the linear scale of the HEALPix 
resolution 512 pixels, is shown in red. Note that the points in 
black have had bandpowers combined together for higher signal- 
to-noise presentation, the best fit is determined from the non- 
combined angular power spectrum. 



dening mask cut. However, closer to the Galactic plane is 
also higher in stripe longitude lambda, so it is possible that 
variable sky brightness could influence the observed galaxy 
density. 

In order to identify the issue causing this underden- 
sity of late-type galaxies, we need to first determine how 
to distinguish these various, nearly degenerate, possible 
causes. First, we extended our algorithm to calculate a cross- 
correlation angular power spectrum. By using this new code, 
we computed the cross power spectrum between late-type 
galaxies and stars, early-type galaxies and stars, and late- 
type galaxies and reddening. From these measurements, we 
find that only the late- type/star cross power spectrum shows 
significant anti-correlation at small t. This evidence implies 
high stellar densities are associated with the slight drop in 
density of late-type galaxies, and suggests that either stellar 
obscuration or star-galaxy separation could be responsible. 

To test the possibility that stars are interfering with 
late-type galaxy densities, we pixelated the stars in the SDSS 
data in the same manner as galaxies and used these to apply 
high stellar density masks to our late- type galaxy sample. 
We found that stringent cuts on high stellar density pixels 
partially reduced, but could not eliminate, this large scale 
power, as shown in black in Figure [TP 

We also attempted to replicate this effect in the early- 
type galaxies, which do not show this large scale power. We 
used the above stellar density masks and artificially low- 
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Figure 12. The ratio of the late-type galaxy radius distribu- 
tion to the early-type galaxy radius distribution as a function of 
Galactic latitude, b. We see that the distribution of galaxy radii 
change significantly at low latitudes. 



Figure 11. The angular power spectrum of late-type galaxies, 
in blue, compared with the same sample masked for high stellar 
overdensities in black. We see that the power in the first band is 
significantly reduced, but still higher than expected. Also in this 
plot, the angular power spectrum of early-type galaxies, in red, 
compared the the same sample with densities artificially lowered 
in high stellar density pixels in green. This demonstrates that 
anti-correlation in star and galaxy densities can result in this 
large scale power, leaving the rest of the angular power spectrum 
largely unaffected. 



ered the early-type galaxy density in those pixels that have 
a higher density of stars. Due to the increased density of 
stars near the Galactic plane, this affects nearly all pixels 
at low Galactic latitudes and high lambda. However, we see 
in green in Figure [TT] that this does raise the large scale 
power in the first bandpower in the same manner as the 
late-type galaxy samples, suggesting a correlation between 
the increased density of stars and the lower than expected 
late-type galaxy density. 

With these results in mind, we next investigated how 
the properties of the late- type galaxy distribution changes as 
a function of Galactic latitude. By comparing the measured 
size of SDSS galaxies u sing the Petrosia n radius containing 
half the Petrosian flux fPetrosian f976), we generated nor- 
malized 2D histograms of of the late- and early-type galaxy 
radii as a function of Galactic latitude. We plot the ratio 
of these normalized distributions in Figure 1121 which shows 
that the galaxy radius distribution changes at low Galactic 
latitudes. The late-type galaxies at low latitudes are on av- 
erage larger than at high latitudes, which when combined 
with the overall reduction in low latitude late-type galaxy 
density suggests that our samples have undercounted small 
late-type galaxies in high stellar density fields. Our inter- 
pretation of these results is that small late-type galaxies in 
the SDSS are being misclassified as stars due to deblending 
issues in crowded fields. 

To ensure that this systematic does not affect our \ 2 



fits for cosmological parameters and bias, we have not in- 
cluded the first bandpower of the angular power spectra de- 
rived from our volume-limited sample in the fitting process. 
We have excluded the first bandpower for all volume- limited 
sample fits for consistency. 



7.2 Comparison with Previous Results 

The constraints provided by our fits to nonlinear theoretical 
matter power spectra for our large multipole and volume- 
limited sample and subsamples are much improved over the 
linear fits used in HBR12. We find that our results are gen- 
erally consistent with an fl m ~ 0.27 with errors on the or- 
der of 0.03, which is typical of galaxy angular power spectra 
dBlake et aZ.ll2007h. This ag rees well with the WMAP7 CMB 
results ijLarson et a/jEmf ) of D. m = 0.267 ± 0.026 as well 



as the combination of SDSS DR8 luminous galaxy angular 
power spectrum resul ts with WMAP 7 and supernova data 
fl m = 0.267 ± 0.0163 i|Ho et aZ.ll20T2l ). 

We find Oj, w 0.03 with errors generally about 0.02 
in our volume-limited and large multipole samples, which 
is consistent with the Q m = 0.0449 ± 0.0028 constraints 
produced by WMAP7. The errors on our measurements are 
an order of magnitude larger than WMAP7, however, as our 
samples are less sensitive to this parameter, largely due to 
uncertainties in the photometric redshifts, an d this is similar 
in oth er galaxy angular power spectra results (jThomas et al\ 
120101) . 

Measurements of the bias parameter vary with the sam- 
ple under consideration. As we've shown in Section[5TT] bias 
is strongly type dependent, so the ratio of early- and late- 
type galaxies in a sample has a profound effect on the bias. 
Even the relative bias between early- and late-type galax- 
ies has wide variatio ns from a r elative bias of 1. 2 ±0.15 
l|Willmer et qZ.lll998f ) to ~ 1.75 (|Ross et al\ 120061 ). and is 
partially dependent on the cut used to separate the differ- 
ent galaxy types. 

With our assumptions of a flat Universe, and various 
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properties of the initial angular power spectrum and reion- 
ization, these results imply that the mass-energy content in 
the Universe is not dominated by mass, but by another form 
of energy, believed to be dark energy. Even the mass in the 
Universe is not primarily normal, baryonic matter, but col- 
lisionless dark matter. The galaxies that we observe make 
up just a few percent of the mass-energy, but these galaxies 
trace the underlying dark matter distribution with a type- 
dependent bias describing how clustered galaxies are com- 
pared to the underlying dark matter. It is worth noting that 
these implications are quite consistent with the WMAP7 
results, despite relying on an entirely independent measure- 
ment process and data set during a completely different cos- 
mic epoch. 

7.3 Future Work 

This technique can be easily applied to other surveys, such 
as the Dark Energy Survey (DES). The DES is scheduled 
to begin taking data by the end of 2012, but even before 
that data is available, simulations of the DES data can be 
analyzed with our angular power spectrum estimation code 
similar to what was done in T02 with the SDSS Early Data 
Release. To support this and similar efforts, we have made 
our parallelized estimation code publicly availablfl 

We also can see the evidence of baryon acoustic oscil- 
lations in the angular power spectrum for narrow redshift 
slices, as in the theoretical angular power spectra in Figures 
[3] and[8] Though the variations in the angular power spec- 
tra from BAOs is generally smaller than our error bars for 
these samples, detection of the BAOs in SDSS galaxy angu- 
lar power spectr a in the 30 < I < 300 range is an area of 
current research (|Seo et aZ.ll2012l ). In deed, the Baryon Qscil - 
lation Spectroscopic Survey (BOSS: fEisenstein et <aZ.|[201lj ) 
of SDSS-III aims to explore the BAO signal through exam- 
ining correlations of Luminous Red Galaxies (LRGs) and 
quasars. 

A closely related measurement that can be done is the 
cross-power spectrum estimation. The method discussed in 
Section [221 nas referred only to auto-correlation, that is cor- 
relating galaxy overdensity with itself; but Equation [2] can 
be extended to calculate the cross-correlation, for example 
between quasars and LRGs, by taking the outer product of 
two different data sets rather than of the same data set. This 
will not only provide interesting scie nce, such as const raints 
on local primordial non-Gaussianity (jSlosar et aZ.ll2008r i. but 
can also be a test of systematics by calculating the cross- 
correlation of two samples that should be uncorrelated, as 
we have seen in Section |7. II 

The advances in general purpose graphics processors 
and their incorporation into modern supercomputers sug- 
gests it is time to revisit the possibility of accelerating this 
computation on these machines. The computationally inten- 
sive parts of this calculation are the signal matrix construc- 
tion and the matrix multiplications and inversions required 
by the quadratic estimator, and we have previously found 
that the se operations perf orm quite well on graphics pro- 
cessors jHaves et al\ I2007T ) . By accelerating these calcula- 

2 We have made these codes freely available at 
http:// ledm. astro. Illinois .edu/ code / apscode . html 



tions, it might be possible to improve the resolution while 
not sacrificing survey area, permitting angular power spec- 
trum measurements with both high bandpower resolution 
and results extending to smaller scales. This could improve 
cosmological parameters constraints from these galaxy sur- 
veys. 



8 CONCLUSIONS 

We have used a quadratic estimator with KL-compression 
to determine the angular power spectrum of the SDSS DR7. 
We applied this technique to a ~ 53.7 square degree area 
near the North Galactic Pole to determine the angular power 
spectrum of the SDSS DR7 main galaxy sample out to high 
multipoles (£ = 800) . By utilizing CAMB to produce nonlin- 
ear power spectra, we created theoretical nonlinear angular 
power spectra by projecting these CAMB spectra to two di- 
mensions using the photometric redshifts of the SDSS data. 
By minimizing the \ 2 fits between these theoretical angular 
power spectra and our measured SDSS angular power spec- 
tra, we find that the best-fit cosmological parameters to this 
SDSS DR7 main galaxy sample are Q b = 0.045 ± 0.012 and 
tt m = 0.267 ± 0.038, with a linear bias of 6 = 1.075 ± 0.056. 

Using the same process on a volume-limited sample cov- 
ering the Northern contiguous SDSS footprint out to £ = 
100, we determine Q b = 0.041 ± 0.020, fi m = 0.282 ± 0.026, 
and b = 1.545 ± 0.057, while the linear bias varies from 
b = 1.515 ± 0.057 in galaxies from 0.3 < z < 0.4 to 
6 = 1.421 ± 0.054 in the range 0.0 < z < 0.3. Splitting 
the volume-limited sample by type, we calculate the relative 
bias b e /bi = 1.375 ± 0.076, and combining this information 
with the photometric redshift slices, we see that the relative 
bias varies from b e /bi = 1.200 ± 0.071 in the 0.3 < z < 0.4 
range to b e /k = 1.421 ± 0.083 when 0.0 < z < 0.3. 

Finally, our examination of the angular power spectrum 
by galaxy type has shown unexpected large scale power in 
late-type galaxies. This suggests that late-type galaxies near 
the Galactic plane are undercounted, and this is associated 
with areas of high stellar density. Our analysis accounted for 
this by excluding the largest scales from the fits to determine 
cosmological parameters. 
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